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Abstract 

We discuss the application of a theory of large-amplitude collective motion 
to a simple model mimicking the pairing-plus-quadrupole model of nuclear 
physics. 

I. INTRODUCTION 

In order to describe processes in nuclei involving large excursions from equilibrium, such 
as shape-coexistence or fission, we need to go beyond the simple harmonic oscillator picture 
underlying the random-phase approximation (RPA). One method is to use a method that 
self-consistently selects collective coordinates for large amplitude collective motion (LACM). 
Recently we have analysed properties of collective motion in a pairing model Hamiltonian 
with level crossings [|I|, where we used a version of a theory proposed by Klein et al 0. to 
determine the collective coordinates. We have shown that the system automatically selects 
either diabatic or adiabatic collective surfaces according to the strength of the pairing inter- 
action. The model used in that paper was rather artificial, since it involved a macroscopic 
system coupling to a few fermionic degrees of freedom. We felt it desirable to test the theory 
for a fully microscopic Hamiltonian which is able to describe systems from vibrational nuclei 
to well-deformed nuclei. 

In this paper, we investigate the collective motion in a microscopic model which describes 
a system of nucleons interacting through a simplified version of the pairing plus quadruple 
force 1^. Although the Hamiltonian has a very simple form, we shall see that the model can 
represent several of the interesting phenomena observed in real nuclei. 

II. THEORY 

The theory of adiabatic large amplitude collective motion (ALACM, see reference for 
a complete description) is applicable to a classical Hamiltonian system which has kinetic 
terms only quadratic in momentum. We thus have to start with a truncated Hamiltonian 

tt) = ^B^^n^np + V{0 , a,P = l,---,n, (1) 
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where the mass tensor B"^, in general, depends on the coordinates and is defined by 
truncation of the Hamiltonian to second order. 

Collective coordinates and intrinsic (non-collective) coordinates g° which are approx- 
imately decoupled from each other, are assumed to be obtainable by making a point trans- 
formation, 

q' = f{0 {^ = h■■■,K), q'' = nO {a = K + 1, ■ ■ ■ ,n) . (2) 

Requiring approximate decoupling between the collective and non-collective coordinates, we 
find a set of coupled equations. In the case of a single collective coordinate {K = 1), the 
equations can be written as 

= A/i, B^W,^,f^^ = uj'fl . (3) 

The equations can be solved iteratively, starting from a stationary point. In principle, 
the procedure to find a collective path is to find successive points at which an eigenvector 
of the covariant RPA equation (second equation) satisfies the force condition (first equation) 
at the same time. 



III. THE MODEL 

The model Hamiltonian we discuss in this paper has 0(4) symmetry.lt has been orig- 
inally developed to describe K'" = 0+ excitations in deformed nuclei 0. It has also been 
used to test various methods for calculation of collective excitations [^-0. The generalisa- 
tion to multiple shells has also been formulated, in order to investigate shape-coexistence 
phenomena 

Consider a single shell with angular momentum j, and define the four basic operators, 

m>0 r)i>0 

^ = X/ ^In^rn , Q = X/ '^mClnCm ■ (5) 
m m 

Here am = —1 if l^^l < ^/2, and +1 for the other values of m, which mimics properties of 
the matrix elements of the quadrupole operator r'^Y^. Q = j + 1/2 is assumed to be an even 
integer. Hereafter we call Q the quadrupole operator. The four operators, P, P, N and Q, 
close under commutation, and generate the Lie algebra so(4). 

We now build a model Hamiltonian with 0(4) dynamical symmetry 

H = -GP^P - \xQ' , (6) 

which contains both a monopole pairing and a quadrupole-quadrupole (QQ) interaction. 
The exact solution can be obtained by diagonalisation in an 0(4) basis. 

The behaviour of the model is based in the competition between the pairing force which 
tends to align quasi-spins and the QQ force which tends to de-align them. This can be made 
extremely clear in the limit of no pairing or no QQ force, where the model can be solved 
exactly. 
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FIG. 1. Excitation energies and transition matrix elements \{n\Q\n — 1)\ (numbers next to 
arrows) in the two-shell 0(4) model discussed in the text, as functions of quadrupole force strength 
X- The case (a) represents the almost harmonic motion for % = 0.1, (b) the flat potential for % = 0.3 
and (c) shape coexistence for % = 0.4. In each case we compare an exact diagonalisation with the 
result of requantisation in a single collective variable. 



We can easily generalise the previous model into a multi-j-shell case, {ji,j2, ■ ■ ■ , Ja)- For 
each j-shell, we take fij = jj + 1/2 to be even, and we weigh the quadrupole operators by qj, 
the magnitude of quadrupole moment carried by the single-particle state j. We again use 
the pairing plus quadrupole Hamiltonian, but we now add spherical single-particle energies, 

H^Y: ^J^mCjm - GP^P - \xQ' . (7) 
jm 



IV. RESULTS AND CONCLUSION 

One of the representative calculations we performed is for a two-shell problem. Both 
shells have Q, — 10, but we take gi = 3, ^2 = 1 and ei = 0, £2 = 10. For strong "deformation" 
we expect that the shells mix strongly, and for weak QQ force we expect a simple spectrum. 
In this case, as can be seen in figure 1, we find that as we move from a harmonic spectrum to 
a more complex mixed situation, that we still keep reasonable agreement for both spectrum 
and transition strengths. The values of the QQ strength used were such that the middle case 
exactly captures the turn-over from single to double well. In all these cases one collective 
coordinate appears to be enough. There are some states in the double-well structure that 
are not well described by the current method, which is currently under investigation. 

In conclusion we have shown that our method can deal with both anharmonic vibrations 
and shape coexistence. In future we shall extend this research to the proton- neutron 0(4) 
model, and to the real pairing-plus-quadrupole model. 

This work is supported by a research grant (GR/L22331) from the EPSRC. 
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